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Abstract. Mass loss on the Asymptotic Giant Branch provides the origin of planetary nebulae. 
This paper reviews several relevant aspects of AGB evolution: pulsation properties, mass loss 
formalisms and time variable mass loss, evidence for asymmetries on the AGB, binarity, ISM 
interaction, and mass loss at low metallicity. There is growing evidence that mass loss on the 
AGB is already asymmetric, but with spherically symmetric velocity fields. The origin of the 
rings may be in pulsational instabilities causing mass-loss variations on time scales of centuries. 



1. Introduction 

The Asymptotic Giant Branch (AGB) contains two separate phases. During the early 
AGB, stars experience quiescent shell burning around the inert but growing carbon- 
oxygen core. Later, quiescent hydrogen burning is punctuated by run-away helium flashes 
(Herwig 2005) followed by short phases of (lower-luminosity) quiescent helium burning. 
The helium flashes are called thermal pulses (TPs). During the TP- AGB phase, the star 
begins to experience strong pulsations, and increasing mass loss. Evolution proceeds from 
semiregular (SR) to Mira and eventually to OH/IR star, with mass-loss rates simulta- 
neously increasing from typically 10 _6 M©yr _1 for Miras, to lO -4 M0yr -1 for OH/IR 
stars. The latter show strong self-obscuration by the circumstellar envelope. Carbon-rich 
stars show a parallel sequence but with lower V-band amplitudes, and higher circumstel- 
lar extinction for the same mass-loss rates. 

The mass loss reaches values much higher than the nuclear burning rate. The hydrogen 
envelope therefore diminishes much faster than that the core grows. Once the envelope 
mass falls below a critical value (~ 0.02 M©), the photosphere begins to collapse, mass- 
loss rates reduce and the star enters the post-AGB phase. In this way, the mass loss 
causes a very sudden end to stellar evolution (e.g. Willson 2000). 

In this paper, aspects of AGB evolution related to the mass loss will be reviewed. 

2. Pulsations 

Pulsation properties evolve strongly during the TP- AGB. The semiregulars have peri- 
ods of typically 50-150 days, small amplitudes and (as the name implies) have irregular 
light curves. Mira variables have periods of 150-500 days, amplitudes in excess of 2.5 mag 
in V, and stable light curves. OH/IR stars show periods of 500-2000 days. The amplitude 
is largest at optical wavelengths. The stellar energy distribution peaks in the infrared, 
and the bolometric amplitude is much less than in the V-band, rarely more than one 
magnitude. The large optical amplitude is caused by variable molecular opacity, mainly 
due to TiO and VO bands. The classification can show imperfections. OH/IR stars tend 
not to be observable in the V-band, and are often classified as SR because the ampli- 
tude is only measured at K. Stars with higher C/O ratio show reduced abundances of 
oxygen-rich molecules, and can therefore also show lower amplitudes. 
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The pulsation equation is given by: 



logP = L9491ogii-0.91ogM-2.07 fundamental mode Wood 1990) 
logP = 1.51ogi? - 0.51ogM + logQ first overtone (Fox & Wood 1982), 

(2.1) 

where in the latter case the pulsation constant Q tzs 0.04; the period P is in days and the 
radius R and mass M are in solar units. Mira variables are fundamental-mode pulsators, 
whilst SR variables show overtone pulsations. Fundamental- mode models including tur- 
bulent viscosity (Olivier & Wood 2005) reproduce Mira light curves fairly well. 

A well-defined Mira period-luminosity relation was found from stars in the LMC. The 
Macho, Ogle and Sirius surveys have shown that this relation is only one of several parallel 
sequences. Magellanic Cloud stars show five sequences (e.g. Kiss & Bedding 2003, Ita et 
al. 2004), two of which are identified as overtones, one is the Mira fundamental mode, one 
is an RGB sequence tracing semi-detached binaries, and the last sequence shows periods 
much longer than the fundamental mode which are still unexplained (Wood et al. 2004). 

Long-term monitoring of Mira variables, mostly by amateur astronomers, has given 
light curves covering, in some cases, several centuries. These extended light curves have 
shown that most Miras are stable, but a number of stars show evolving periods. The best 
known case is that of R Hya, for which historical records (Zijlstra et al. 2020 indicate 
a continuous decline from about 495 days before 1800, to 385 days in 1950. The period 
has remained stable since 1950, and the earlier longer period also may have been stable. 



Templeton et al. (2005) has analyzed the full set of AAVSO records available, and out of 
547 Mira variables find 8 cases of long-term period changes, at Qsigma significance. The 
most dramatic case is T UMi, which showed a constant period for over 60 years, before 
a sudden, sharp decrease in period by 25 per cent over the next 20 years. 



Zijlstra & Bedding (2002) define three types of period instability: sudden changes as 



in T UMi, continuous changes as in R Hya, and meandering periods. The latter appears 
to be common among the longest period Mira variable, and involve variations up to 10 
per cent on time scales of ~ 50 years. 

Several explanations exist for the period changes. Thermal pulses involve large, sudden 
variations in luminosity (and therefore radius and period). This fits stars such as T UMi 
well, where there was no pre-indicator of a developing instability. The rate of thermal 
pulses is of the order of 0.5 per cent per century, or about ~ 2 per century among the 
known Miras. The second possiblity is a non- linear pulsation (Olivier & Wood 2005), 
where the change in amplitude causes a change in period (Ya'ari & Tuchman 1996). 
Following a large amplitude change, e.g. the onset of pulsation, they find a large period 
decrease of the fundamental mode over a time scale of centuries. The observed correlation 
between period and amplitude changes (e.g Zijlstra & Bedding 2002) provides some 
support for this model. Icket et al. (1992)] have pointed out that the surfaces of long- 



period stars are susceptible to weak chaos. |Zijlstra et al. (2004) discuss strongly variable 
periods among SC stars, where the C/O ratio is very close to unity and the chemical 
equilibrium is highly sensitive to temperature changes: they suggest that in SC stars, 
rapid changes in molecular abundances may trigger weakly chaotic period changes. 



3. Mass loss 

3.1. Proposed relations 

AGB mass loss is a multi-step process. The pulsations extend the atmosphere, and allow 
molecules to form. At large distances the temperature drops below the condensation tern- 
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peratures and a number of solids form: MgSiO-bearing grains in O-rich environments, and 
SiC and amorphous carbonaceous particles around carbon-rich stars. Radiation pressure 
pushes the grains out, and by collisions carries the gas with it. Pulsation can support 
mass-loss rates ~ 10~ 7 M Q yr -1 , but the radiation pressure on dust (and molecules) 
greatly increases the mass loss. The difference corresponds to the case A and case B 



mass-loss regimes of Winters et al. (2000) 



(i) 


logM 


= aP + b 


Vassiliadis & Woodl993 


(ii) 


M 


= aMr 2 - 1 L 31 RM- 1 


Bloecker 1995 


(ii) 


M 


= L2.47 T -6.8 M -1.96 


Wachter et al. 2002 


(iv) 


M 


= a £l.06 T -6.3 


van Loon et al. 2005. 



Dynamical models of the extended envelopes do not yet convincingly predict mass-loss 
rates. Instead, formulistic mass-loss rates are used. The following relations are in general 
use: 



(3.1) 



Relation (ii) is derived by fitting the Bowen pulsation models to observed initial-final 
mass relations; (iii) is derived for carbon-rich stars only, using dynamical models, whilst 
(iv) comes from mass- loss determinations in the LMC, including both AGB and RSG 
stars, and is observationally the best established. 

A very strong dependence on stellar temperature is indicated. However, the optical 
spectral type used to assign a temperature measures the r = 1 layer in the molecular 
envelope, far above the real stellar surface: it underestimates both the surface and the ef- 
fective temperature. The strong dependence is also deceptive. Using L oc R 2 T 4 , relations 
(iii) and (iv) give M oc T 36 and T -2 - 3 , respectively. The first one is physically unlikely, 
and Wachter et al. (2002) may overestimate the effect of luminosity. The Bloecker rela- 



tion suffers from the same problem. Relation (iv) excludes mass as a relevant parameter, 
because this parameter is not known for the LMC stars. However, we may adopt the 
suggestion that the mass loss depends on the surface binding energy, M/R. In this case, 
relations (iii) and (iv) suggest that the total dependence is along the lines of 

M ~ (j-T) IP, (3.2) 

with a ~ —2; the luminosity dependence is likely weaker, with 7 =~ 1. 

In addition to the surface parameters, the mass loss is also affected by the period, 
since the atmosphere extension is pulsation-driven. Relation (i) is the only one to include 
the pulsation explicitly. But the other relations also implicitly do so, via the pulsation 
equation. However, they predict a period dependence which is much weaker than found 
for (i). An example to constrain the period dependence is given by the evolution of R 
Hya, because its decreasing period over the past 200 years was accompanied by a strong 
decrease in mass loss (Zijlstra et al. 2002). Its evolution fits relation (i) well, but is not 
be reproduced by the other relations. 

Overall, the predictive value of the proposed relations is still limited. Important pa- 
rameters, such as composition and metallicity, are not or not fully included. The case A 
and B regimes may show different dependences on stellar parameters. 

3.2. Time-variable mass loss 

Mass loss is known to vary on several time scales. Over 10 6 yr, the slow increase in the 
absolute magnitude of ~ 1 mag will be accompanied by an increase of M e.g. Vassiliadis 
& Wood 1993. This slow evolution is interrupted by the thermal pulses, which cause a 
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sharply enhanced mass loss event, occuring every 10 to 10 5 yr, depending on the core 
mass. Mass loss variations on a time scale of 10 2 -10 3 yr are implied by the multiple rings 
seen around many PNe. Finally, extinction variations in both carbon-rich and oxygen-rich 
stars indicate changes within decades. 

Mass loss spikes during the TP are the accepted explanation for the detached shells 
seen around TT Cyg, U Cam and DR Ser, among others (e.g. Olofsson et al. 2000); the 
thin shells can be explained by a brief burst of 100-fold stronger mass loss, sweeping 
up the previous, weaker and slower wind (Schoier et al. 2005). The TP is followed by a 
subsequent phase of quiescent helium burning when L is a factor of 3 lower than during 
the hydrogen burning phase, with a weaker wind. 

The multiple, thin rings seen around PNe (Corradi et al. 2004) still lack a commonly 
accepted explanation. First, an instability involving the dust driving was suggested by 



Simis et al. (2001) dust formation leads to increased dust density, until dust-gas cou- 
pling empties the dust shell. This leads to mass-loss episodes at the correct time scales. 
However, this instability is not correlated between different directions, and does not lead 
to spherical structures. Cooling due to existing dust strongly enhances the local dust for- 
mation rate, leading to near-chaotic mass loss Woitke (2005) Second, a nuclear-burning 



instability was sugg ested by |Van Hron et al. (2003)| Third, the period instability found 



for long-period Miras, interpreted as due to a non-linear pulsation, could cause the rings 
(Zijlstra & Bedding 2002, via period-induced mass-loss variations. The hypotheses need 
to explain the slight irregularity in the ring spacings, and the decrease in spacings for 
rings closer to the star (i.e. ejected later). The model involving the non-linearity in the 
pulsations can explain both. However, the issue is still far from settled. 

The decadal (extinction) variations are most easily explained by the dust clumping 
as found in the Woitke models. Extinction variations are mainly seen in carbon-rich 
stars (e.g. Whitelock et al. 2006) but a few oxygen-rich stars also show variations due to 
circumstellar extinction, most noticeably L2 Pup (Bedding et al. 2002). 



4. Morphology 

Strong asymmetries are present in the morphologies of almost all post-AGB envelopes. 
These structures must have their origin during the AGB mass loss. But remarkably little 
evidence for the precursors of the post-AGB morphologies has been found. In fact, the 
current consensus is that AGB envelopes are spherically symmetric, and that asymmetries 
form during a run-away process at the very tip of the AGB. The argument that such a 
fundamental change happens during a phase which is unobservable only by its brevity, 
is far from convincing. 

Two observations provide evidence for spherical AGB shells. One comes from the OH 
maser profiles of OH/IR stars, with strong peaks at the extreme velocities. For asym- 
metric envelopes, the emission at intermediate velocities becomes stronger and multiple 
peaks may develop (Chapman 1988). The second observation is that the outer structures 
around PNe and post-AGB stars are often round, especially so for the rings. The regular 
detached shells around, e.g., TT Cyg should also be mentioned. 

However, both can be questioned. Maser emission cannot be used to determine a 
density structure. Instead it is an exponential function of the column density at constant 
velocity, in competition between all directions. OH profiles are far more dependent on 
the velocity structure than on the densities. An asymmetric envelope with a spherically 
symmetric velocity structure will give rise to symmetric profiles. 

The rings appear to provide stronger evidence. However, any mass-loss changes leading 
to such rings are likely to be accompanied by changes in the expansion velocity. The 
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Figure 1. NTT V-band image of RZ Sgr, showing its bipolar reflection nebula 
rings then become the swept-up interface between the more massive, faster wind and 



the slower, weaker wind: Zijlstra et al. (2001) show that the swept-up shell will move at 



constant velocity, with the location at anyone time given by: 



^ring^ — ^slow 



l-M+(£-l)V/Vl , (4-1) 



where 



''"fast t Vfast 

a* = ; £ 



If the two winds are self similar, and specifically if they show spherically symmetric 
velocities and show the same ratio of mass-loss rates in all directions, v r i ng will be direction 
independent. The swept-up interface will be spherical, even if the underlying mass loss 
is not. 

Thus, in both cases there is strong evidence for spherically symmetmcvelocity fields, 
but the mass loss itself can be asymmetric. 

There is in fact growing evidence for bipolar shells on the AGB. The bipolar enve- 
lope of X Her was mapped by |Nakashima (2005)] The reflection nebula around RZ Sgr 
Whitelock (1994]] is bipolar in our NTT image (Fig. ^| . Linear polarization is present in 



integrated photometric data of some OH/IR stars (Jones & Gehrz 1990). Finally, the core 
of IRC+10 216 shows several components, possibly in a bipolar configuration (Osterbart 
et al. 2000). (The last two points can also be expained by Woitke-type blobs (Woitke 
2006), however.) 



5. Binarity 

Whereas binarity is now thought to be very common among post- AGB stars and PNe, 
and the likely cause of the asymmetries of the envelopes, very little is known about 
binarity among AGB stars. The extreme brightness and low temperature of the AGB 
star make a binary comnpanion difficult to detect, except where the companion interacts 
with the extended atmosphere or the wind. Direct evidence for binarity among AGB stars 
will be needed to make evolutionary connections between AGB and post-AGB stars. 

A distinction can be made between four different types of interaction. 

(a) Common envelope systems, where the companion triggers the mass loss. These 
system will rarely reach the AGB, as the common envelope is likely to develop already 
on the RGB. Orbital separations are of the order of 1 AU. Perhaps 10 per cent of 'proto 
planetary nebulae' may derive from such systems, although fewer would evolve into PNe. 
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(b) Wider systems, where the geometry of the mass loss is affected, but not the mass- 
loss rate itself. These systems may include symbiotic stars and the RV Tau progenitors. 
Orbital separations are a few to tens of AU, and the systems exhibit circumbinary disks 
and accretion. These may account for ~ 25 per cent of the PN birth rate. The companion 
may ionize part of the stellar wind, as in OH 231.8+4.2. 

(c) Systems with only minor effects on the mass loss geometry. Typical separations 
are ~ 100 AU. The geometry may show the movement of the mass-losing star, as in the 
spiral nebula around AFGL 3068 (Mauron & Huggins 2006, see also Morris et al, these 
proceedings). Mira itself has a white dwarf companion ~ 100 AU away. 



6. ISM interaction 

The study of the interaction between the AGB wind and its local ISM was pioneered 
by Villaver et al. (2002) Previously, the interaction had only been considered for old, 
ionized PNe. The interaction begins in fact with the onset of the AGB wind, and leads 
to a swept-up shell surrounding the AGB star. Once the swept-up mass exceeds the 
mass lost in the wind, the shell slows down and eventually becomes near-stationary. The 
typical radius of this shell is 1-2 pc, depending on local conditions. An example of such 
a shell, called a 'wall', was found by Zijlstra & Weinberger (2002) The wall may wipe 
out any memory of long-term mass-loss fluctuations. Once the PN fast wind reaches the 
outer wall, the PN rebrightens; this increases the PN life time. 

If there is a significant proper motion of the star with respect to the ISM, the wall 
becomes one-sided and at higher velocities, becomes a bow shock. This is modelled by 



Wareing et al. (2006) based on the IPHAS morphology of the PN Sh2-188. The Spitzer 



image of R Hya (Speck, these proceedings) show another example of such a bow shock.) 
Evolutionary calculations for a large range of parameters are presented by Wareing (these 
proceedings). 



7. Metallicity 

The effect of metallicity on AGB mass loss is only just beginning to be explored. 



Bowen & Willson (1991) first showed that for [Fe/H]< —1, dust-driven winds fail and 



the mass loss becomes pulsation-driven instead. This significantly reduces the mass- 
loss rates. Using their formalisms, Zijlstra (2004) find that the final (white dwarf) mass 
becomes higher at very low metallicity, and may reach the Chandresekhar limit. This 
does not happen in the models of Herwig (2004) where core-dredge-up limits the core 
mass to k 1 Mq . 



Mass-loss rates in the LMC are as high as those in the Galaxy Van Loon et al. (1999) 



However, it is possible that these values are reached at higher luminosities than for Galac- 
tic stars. (The metallicity of the LMC is comparable to the outer Galactic disk.) Evidence 
for reduced mass-loss rates at much lower metallicity comes mainly from the specific fre- 
quency of PNe, which is significantly lower in metal-poor dwarf galaxies compared to 
more metal-rich systems. This is illustrated in fig. [21 

Expansion velocities will be much lower at lower metallicity, limited by the dust-to- 
gas ratio. Evidence is sparse, as expansion velocities outside of our Galaxy have only 
been measured for some LMC OH/IR stars. These indeed indicate a significantly lower 
expansion velocity in the LMC (Marshall et al. 2004). 

Dredge- up of carbon will lead to higher C/O ratios at lower metallicity, because of the 
lower oxygen abundance. This leads to the counter-intuitive effect that carbon molecules 
are more abundant at lower metallicity. The C2H2 infrared bands were indeed found to be 
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Figure 2. Left: The ratio between number of planetary nebulae and luminosity of the 
parent stellar popula tion (i/-axis), as function of metallicity (a;), for nearby galaxies, from 
Magrini et al. (2003) | Encircled points indicate spiral galaxies affected by internal extinction. 
Right: The same ratio, as function of the ratio of carbon over M-type stars on the AGB, a tracer 
of the metallicity of the AGB population 



very strong in carbon stars in the Magellanic Clouds (e.g. Matsuura et al. 2005). This has 
consequences for dust formation: amorphous dust may form easily around low metallicity 
carbon stars, while minerals such SiC and MgS do not, ginvig a mineral-poor or 'soft' 
dust. Sloan et al. (2006) and Zijlstra et al. (2006) using Spitzer spectra, indeed find that 
the SiC and MgS features are weaker in the Magellanic Clouds. Oxygen-rich stars can 
only form dust from metal-dependent minerals, and their dust formation efficiency is 
expected to be strongly suppressed at low metallicity. Mass loss and dust formation at 
low metallicity is therefore expected to be dominated by carbon stars. 
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